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Microquasars are X-ray binaries with relativistic jets. These jets are powerful energy 
carriers, thought to be fed by accretion, which produce non-thermal emission at different 
energy bands. The processes behind the bulk of the non-thermal emission in microquasars 
may be of leptonic (synchrotron and inverse Compton) and hadronic (proton-proton in- 
teractions, photo-meson production, and photo-disintegration) nature. When leptonic, 
the fast particle cooling would allow the obtention of relevant information about the 
properties close to the accelerator, like the radiation and the magnetic field energy den- 
sities, and the acceleration efficiency. When hadronic, the extreme conditions required in 
the emitter would have strong implications on the physics of jets and their surroundings. 
The very high-energy part of the spectrum, i.e. > 100 GeV, is a good energy range to 
explore the physics behind the non-thermal radiation in these compact variable sources. 
In addition, this energy range, when taken altogether with lower energy bands, is a key 
piece to construct a comprehensive picture of the processes occurring in the emitter. 
Until recently, the very high-energy range was hard to probe due to the lack of sensitiv- 
ity and spatial and spectral resolution of previous instrumentation. Nowadays, however, 
powerful gamma-ray instruments are operating and the quality of their observations is 
allowing, for the first time, to start to understand the production of high-energy emission 
in microquasars. 

To date, several Galactic sources showing extended radio emission, among them at 
least one confirmed microquasar, Cygnus X-l, have shown a TcV signal. All of them 
show complex patterns of spectral and temporal behavior. In this work, we discuss the 
physics behind the very high-energy emission in the microquasar Cygnus X-l, and also in 
the other two TeV binaries with detected extended outflows, LS 5039 and LS I +61 303, 
pointing out relevant aspects of the complex phenomena occurring in them. We conclude 
that the TeV emission is likely of leptonic origin, although hadrons cannot be discarded. 
In addition, efficient electromagnetic cascades can hardly develop since even relatively 
low magnetic fields suppress them. Also, the modeling of the radiation from some of the 
detected sources points to them as either extremely efficient accelerators, and/or having 
the TeV emitter at a distance from the compact object of about ~ 10 12 cm. Finally, 
we point out that the role of a massive and hot stellar companion, due to its strong 
photon field and wind, cannot be neglected when trying to understand the behavior of 
microquasars at high and very high energies. The complexity of microquasars precludes 
straightforward generalizations to a whole population, and are better studied presently 
in a source by source base. The new and future gamma-ray instrumentation will imply 
a big step further in our understanding of the processes in microquasars and gamma-ray 
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emitting binaries. 
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1. Introduction 

Microquasars are an X-ray binary (XRB) subclass formed by those sources that 
present extended radio jets (e.g. Mirabel & Rodriguez 0}. These systems are formed 
by a non-degenerated star, which can be in different stages of its evolution, and a 
compact object, which can be a black-hole or a neutron star. Depending on the mass 
of the non-degenerated stellar companion, the system is considered a low or a high- 
mass microquasar. Typically, systems harboring an OB star are considered high- 
mass XRBs (HMXB), and XRBs with later type stellar companions are classified as 
intermediate or low-mass XRBs. It is thought that the compact object powers the 
relativistic jets via accretion of matter expelled from the companion. This material, 
when reaching the surroundings of the compact object, forms an accretion disk that 
is usually detected in the X-rays. Simplifying very much the case, depending on the 
accretion state, the X-ray spectrum varies strongly, from a multi-color black body 
peaking around 1 keV plus a minor soft power-law spectrum at higher energies, 
for high accretion rates (high-soft state), to one dominated by a hard power-law 
spectrum plus an exponential cutoff around 100 keV, for low accretion rates (low- 
hard state; for an extensive description of the X-ray phenomenology, see McClintock 
& Remillard^). It is expected that a persistent jet will be present during the low- 
hard state, and a transient ejection will form when switching from the low-hard to 
the high-soft state (Fender et al.^J). Correlations between the radio and the X-ray 
luminosity, and the accretion/ejection activity, have been proposed (e.g. Gallo et 
al.Sl, Corbel et al.El, Fender et al.% 

The jet formation and the production of non-thermal radiation in the jet are 
major ingredients that distinguish a microquasar among other types of XRBs. The 
non-thermal radiation produced in microquasar jets has been resolved in radio at 
very different spatial scales (e.g. Ribo ^) and also in X-rays at large scales (e.g. 
Corbel et al. ®). This emission is a clear evidence that particle acceleration takes 
place in different locations of microquasar jiets under very different conditions (for 
a discussion on this, see, e.g., Bosch- Ramon^j). There are also radio emitting XRBs 
in which extended emission has not been detected. I t ha s been proposed that these 
sources could be microquasars as well (e.g. Fender 0^2), like the low-mass system 
XTE J1118+480 (e.g., Chaty et al.Hj), the jet of which has not been resolved yet. 

Microquasars show up themselves as compact and rapidly variable sources from 
radio to very high energies. In such a type of emitters, when radiating particles 
are leptons, the highest energy part of the spectrum is a good range to explore 
non-thermal processes. It is due to the short timescales associated to the particles 
that produce the emission, which implies that the accelerator and the emitter are 
likely the same or similar regions. In case of hadrons, it is possible to derive impor- 
tant information concerning the jet hadronic content, at least of its relativistic part, 
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whereas at the same time it is giving information of the conditions of the emitter, 
like very dense matter and/or target photon fields. Moreover, photons generated 
by very high-energy (VHE) electrons and/or protons give a better insight on the 
mechanism of acceleration and the conditions under which it takes place, helping 
to understand better the processes that accelerate particles up to such high ener- 
gies. Finally, the presence of a hot and massive star renders a scenario in which 
photon-photon absorption, and the occurrence of electromagnetic cascades, can be 
studied. This can give important information on the conditions of the massive star 
surroundings. 

Historically, the poor spatial resolution and sensitivity of the available instru- 
ments working at gamma-ray energies were not enough for accurate theoretical 
modeling, although these sources were proposed to be gamma-ray emitters more 
than a decade ago (e.g. Levinson & Blandford^^l Levinson & MattoxQ21 Paredes 
et al. ^j). Nowadays, however, powerful gamma-ray instruments are operating or 
will start soon, and the quality of their observations is allowing us, for the first 
time, to probe the physical processes that take place in microquasars and their jets. 
These observations at very high energies give the necessary input to constrain the 
theoretical models that lacked in the past. 

The recent evidence of detection of a transient ev ent, at the (post-trial) signif- 
icance level of 4.1 a, from Cygnus X-l (Albert et al. have shown that micro- 
quasars can indeed produce VHE emission, and VHE observations can bring to us 
important information on these sources. Another two interesting cases are LS 5039 
and LS I +61 303. The former has been detected in the TeV range by the Cherenkov 
telescope HESS (Aharonian et al. H^j) , showing a periodic behavior in the VHE ra- 
diation (Aharonian et al. with the same period as the orbit (Casares et al. ^^J) . 
LS I +61 303 has been also detected at TeV energies by the Cherenkov telescope 
MAGIC, being its emission variable (Albert et al. 1^). Unlike Cygnus X-l, which is 
a firmly established microquasar, LS 5039 and LS I +61 303 are no considered at 
present microquasars due to some peculiarities in their X-ray and radio characteris- 
tics. Nevertheless, both sources present, like Cygnus X-l, extended radio emission, 
and harbor OB stars. 

In this review, we want to take profit of the new microquasar phenomenology at 
very high energies, to draw a theoretical picture of the non-thermal processes that 
could take place in these sources. This is to be put in the context of the historical 
evolution of the field, which will be summarized. To explore which are the relevant 
processes in the microquasar scenario, we will carry out a detailed review of different 
mechanisms: particle transport, radiation and photon-photon absorption, in the 
context of microquasars. Using this sound and basic theoretical background, plus 
the observational knowledge at very high energies from the sources presented above, 
the microquasar Cygnus X-l, and also the TeV binaries LS 5039 and LS I +61 303, 
constraints on their physical conditions will be inferred. We note that the type of 
approach applied to these sources is applicable to large extent to any close binary 
system emitting in the TeV regime. It is worth mentioning here a recent review 
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of a broader topic by Levinson (J^) on VHE radiation from jets including active 
galactic nuclei, gamma-ray bursts, and microquasars. It is also interesting to note 
that massive young stellar objects, also presenting jets, have been proposed to be 
gamma-ray emitters by Araudo et al. 

This work is organized in the following manner. In Sections [21 and [31 we try 
to look comprehensively at the microquasar phcnomcnological picture, setting up 
a general scenario for these sources. In addition, there we summarize in a non- 
exhaustive way previous studies of different topics related to microquasars at high 
energies. In Sect.|4j a short review on different particle transport, acceleration and 
radiation mechanisms is done, treating also briefly the issue of gamma-ray absorp- 
tion and electromagnetic cascading in microquasars. All this can help to understand 
how the new VHE data fits in previous ideas and frameworks. Also, we put forward 
a plain but physically sound leptonic model to explore the processes that are in- 
volved in the generation and absorption of VHE emission in jet galactic sources, 
trying to understand which kind of physics is relevant there. For this, the recent 
observational findings at very high energies concerning microquasars are used. In 
Sect. [5J some hot topics of the field are discussed with some extension. Finally, in 
Sect. [51 the main ideas of this work are summarized. 

2. The study of microquasars 

At the discovery of microquasars, their physics was object of speculation due to 
the morphological simil ariti es of these sources with their larger scale analogs, the 
quasars (Mirabel et al. ^J). The question was to find out to which extent the 
processes taking place in extragalactic jet sources could be extended to galactic 
ones. In addition to morphology, kinematical similarities became apparent when 
superluminal movi ng ejecta were found at radio wavelengths in GRS 1915+150 
(Mirabel et al. ED). Nowadays it is widely spread the opinion that galactic and 
extragalactic relativistic jet sources share much more than just morphological and 
kinematical resemblance. In fact, since both types of source harbor a compact object 
surrounded by an accretion disk and a relativistic outflow that is an efficient non- 
thermal emitter, the original analogy argument has been extended not only to the 
mechanisms that are producing such non-thermal radiation but even beyond. The 
aim would be to embrace also the physical link between the relativistic particle 
generation, the accretion phenomena, the jet formation and to some extent the 
interaction between the jet and its environment. A well known example would be 
the phenomenological scaling laws relating the physics of the radio and the X-ray 
emission in microquasars with the accretion rate (e.g. Corbel et al.^j Gallo et al.^j), 
and its more general version that includes t he m ass of the bla ck- h ole, being exte nde d 
to extragalactic objects (e.g. Merloni et al. ^2 Falcke et al. Kording et al. ^S). 
As already mentioned, these empirical scaling laws include the accretion rate and 
the central object mass as the key parameters to account for (at some stage Doppler 
boosting is to be added), as well as phenomenological links, sometimes accretion and 
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jet formation theory motivated, between the radio, the X-ray, and the jet kinetic 
luminosities and the accretion rate itself. 

It is usual and seems natural to associate the properties of the accretion disk, the 
jet and the compact object mass. The latter will influence strongly the properties 
of the infalling matter, like the accretion rate, the accreted matter temperature, the 
magnetic field strength or the accretion disk radiative output. In many XRBs, stan- 
dard accretion theory can explain the X-ray spectrum quite accurately via thermal 
emission from the inner parts of the accretion disk plus a corona-like emitteiQ, al- 
though there is an on-going debate concerning the origin of the hard X-rays, whether 
they come from the jet base, or from a corona-like reg ion close to the compact ob- 
ject, both of similar properties (e.g. Markoff et al . EST Maccarone Efijj j n an y casej 
some sources do not fit in any of these schemas, since their X-ray radiation does 
not seem to come from the regions close to the compact object but further out, like 
the case of SS 433 and, perhaps as well, LS 5039 and LS I +61 303. This could be 
explained by recent theoretical studies that show that the role of the magnetic field 
can be more crucial than the compact object mass itself and can lead to radiatively 
inefficient accretion (Bogovalov & Kelner , in which case the jet could dominate 
in the X-ray band (see also Fender et al. ®) . 

Theoretically, the relationship between the compact object mass and the jet 
physics (e.g. ejected mass rate, matter content, internal energy, carried magnetic 
energy, bulk speed) is unclear since the accretion/ejection physics is not well known. 
Furthermore, the link between the main jet properties and the production of non- 
thermal particles and their radiation cannot be based on first principles, since it 
relies on the particular conditions of the jet plasma (sound speed, diffusion coeffi- 
cients, local magnetic field, degree of turbulence, presence of shocks and their veloc- 
ity, radiation field, plasma density and temperature, etc.), which are not properly 
constrained. Actually, the situation is even more complicated due to the influence 
on the jet processes and the non-thermal radiation itself of external factors like, for 
instance, an external radiation field, the wind produced by the star, the properties 
of the interstellar medium (ISM), etc. Therefore, although empirical laws can help 
to classify objects via significant features of their emission and be powerful heuristic 
tools, it is sometimes difficult to motivate them from first principles, and fundamen- 
tal approaches become necessary. Such fundamental approaches have to be applied 
source by source and grounded on detailed observational data. It requires keeping 
a comprehensive approach to study the source (i.e. phenomenological studies), but 
developing more fundamental theoretical models. 

3. High energy processes in microquasars: a historical perspective 

In this section, a schematic and descriptive picture for the relevant processes taking 
place in microquasars, like particle acceleration, and generation and absorption of 



a A region surrounding the compact object that would be filled with a hot population of electrons 
and ions. 
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VHE radiation, is presented together with a historical perspective of the held. This 
will set up the context within which a more detailed physical treatment will be 
carried out in Sect. HJ 

3.1. A descriptive picture 

For illustrative purposes, we show in Figure [1] a very schematic picture of a micro- 
quasar together with its main elements, and some of the processes taking place in 
these objects. Some regions of the jet are labeled depending on the main radiative 
products. This simplified picture would consist on several important elements: the 
companion star, the stellar radiation held, the stellar wind, the compact object, 
the accretion disk, and the jet itself. The relativistic particles inside the jet will 
interact with the stellar radiation field as well as the jet magnetic field producing 
inverse Compton (IC) and synchrotron radiation, respectively. Under the strong 
stellar radiation field, the creation of pairs will be unavoidable if VHE gamma-rays 
are produced. Finally, the stellar wind may have a significant impact on the jet. All 
this, plus some additional physical processes (not included in Fig. [T] for clarity), like 
hadronic emission, will be discussed in more detail in Sect. [4] 

3.2. Jet formation, evolution and termination 

We present here a brief description of a plausible picture for the jet development, 
from the launching to the termination point. 

Although a complete theory for jet generation is still lacking, many studies on 
jet powering, acceleration and collimation have been carried out during the last 



Fig. 1. A sketch of a microquasar. We explicitly show the star, the stellar wind, the accretion 
disk, the jet, the jet relativistic electrons, some leptonic radiation processes, and pair creation 
under the stellar photon field. 
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decades (e.g. Blandford ^ Blandford & Znajek ^ Blandford & Payne ^ Meier 
|31 Koide et al.^S Chattopadhyay & ChakrabartiESl Meier^ Hujeirat^ Meier 

Bogovalov & Kelner ^ De Villiers et al. I3S1 Ferreira et al. I^McKinney ED 
Hawley & KrolikE21 Komissarov et al. E2 Barkov & Komissarov EH) At present, 
due to the apparent correlation between accretion and jet activity (e.g. Fender et 
al.^, a widely accepted scenario is one in which jets are powered by accretion. The 
accreted matter starts to move following ordered magnetic field lines that thread 
the inner regions of the accretion disk. By magneto-centrifugal forces, the plasma 
is ejected by these magnetic field lines from the accretion disk in the direction 
perpendicular to it. The differential rotation of the accretion disk would create 
a spiral-like shape of the magnetic lines. This magnetic field configuration would 
accelerate and collimate the plasma. Although it is unclear at which scales the 
jet is already formed, VLBI observations of extragalactic jets show evidences that 
the collimation region could be located at ~ 100-1000 i?sch (e.g. Junor et al. E^l 
Horiuchi et al. ESQ 

Once the microquasar jet is formed, it may interact with dense material ejected 
by the accretion disk or the stellar companion, the latter being almost unavoidable 
in high-mass systems. Although the role of an accretion disk wind could be to give 
further collimation and stability to the jet (e.g. Hardee & Hughes SHI Tsinganos et 
al. ESI), the role of a strong lateral stellar wind may lead to jet bending and even 
disruption (Perucho & Bosch- Ramon ^21). In any case, jet /environment interaction 
can lead to shock formation and the radiative counterpart may be observable either 
as a transient phenomena, when the jet (or a discrete ejection) penetrates for the first 
time the surrounding medium, or as a (quasi) steady one, when the jet formation 
is continuous at t he relevant timescales and recollimation shocks occur (Perucho & 
Bosch- Ramon l^2J) . 

Despite the environment of microquasars at large spatial scales may be quite 
different depending on the Galaxy region or the strength of the companion star 
wind, the jets should stop somewhere, terminating either via disruption, in a similar 
way to extragalactic Fanaroff- Riley I sources (e.g. Fanaroff & Riley Perucho & 
Marti 1SZ1) ; or via strong shocks in the ISM, as seems to be the case for Fanaroff- 
Riley II sources (e.g Fanaroff & Riley ^D Kaiser & Alexander ^31 Scheck et al.^H). 
In either case, the radiative outcomes would be different. A classical example of the 
interaction between a microquasar jet and its env ironment at large sc ales is the case 
of the W 50-SS 433 system (e.g. Zealey et al.ESl Velazquez & Raga^. 

In the following sections, we list several non-thermal processes, i.e. particle accel- 
eration and radiation, that may occur in the jet giving a (non exhaustive) overview 
on the literature. For the sake of clarity, we have divided the jet in several regions. 



b The jets of the microquasar SS 433 present collimation di stan ces, inferred from X-ray observa- 
tions, apparently larger, of about ~ 10 6 i?Sch (Namiki et al.H3). We note however that the jet of 
SS 433 is quite different from a typical compact microquasar jet because of its huge kinetic power, 
messy environment, and thermal emitting nature. 
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The jet base, close to the compact object (~ 100 — 1000 Rsch), a region farther out, 
at the binary system scales (typically ~ 10 11 — 10 13 cm), well outside the binary 
system -jet middle scales- (around 10 15 — 10 16 cm), and the termination regions of 
the jet, where it ends interacting somehow with the ISM (> 10 17 cm). We show in 
Fig. [2] a sketch of the different considered regions. Later on, with the most recent 
observational findings at hand, we will go deeper in the modeling to find out which 
of the different considered scenarios is the most suitable to explain the data. 



3.3. Particle acceleration 

Particle acceleration processes in microquasar jets can be of different types. At the 
base, the jet could be magnetically dominated (e.g., in extragalactic jets: Sikora 
et al. , and a mechanism of particle acceleration may be mag neti c energy dis- 
sipation via MHD instabilities in the jet base (e.g. Zenitani et al. EBJ) Also, if jet 
velocities were high enough in the base, the dense available photon (from an accre- 
tion disk/corona) and m atte r fields could allo w th e converter mechanism to take 
place (e.g. Derishev et al. ESI Stern & PoutanenEOl) Magnetic field reconnection in 
the surrounding corona could inject a non-thermal population of particles in the jet 
as well (e.g. Gierlinski & Done EU and references therein). A magneto-centrifugal 
mechanism could also operate very close to a rotating black-hole (e.g. Neronov & 
AharonianE21 Rieger & AharonianESl) 

At binary system scales, plausible mechanisms to generate relativistic particles 
in the jet arc the different versions of the Fermi process: shock diffusive (Fermi I); 
random scattering (Fermi II); and shear acceleration (e.g. DruryEll Fermi ESI and 
Rieger & Duffy ES respectively; see also Rieger et al. E3). Fermi I mechanism 
could take place due to internal shocks in the jet; Fermi II acceleration could take 
place if magnetic turbulence is strong enough, with high Alfven velocity; shear layer 
would be a natural outcome of an expanding jet or different jet/medium velocities. 




10*9 10*11 10*13 10*17 10*18 Cf 

Jet base Binary system scales Jet middle scales Jet termination region 



Fig. 2. Sketch of a microquasar with the different relevant scales resumed in this section. 
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Interactions with the stellar wind may also trigger particle acceleration via, e.g., 
a recollimation shock formed in the jet that expands against the dense material 
expelled by the companion star (e.g. Perucho & Bosch- Ramon . The velocities 
of the shocks mentioned here could be either mildly or strongly relativistic. In the 
latter case, the converter mechanism may be effective in very bright star systems. 

At microquasar jet middle scales, some sort of shock acceleration might still take 
place. For instance, intermittent ejections at timescales of ~hours-days and different 
velocities could create shocks at distances of about ~ 10 15 — 10 16 cm. Also Fermi II 
type and shear acceleration appear plausible for a continuous outflow at these scales 
(something simila r co uld happen in the intra-knot regions of extragalactic jets, see, 
e.g., Rieger et al. ^j). Regarding the environment dynamical influence on the jet, 
it is not expected to be significant given the high jet ram pressure compared with 
the medium one at these scales. 

At the microquasar jet termination point, as in AGN hot spots and radio lobes 
(e.g. Kaiser & Alexander ^2J) , the external medium plays an important dynamical 
role. When the swept ISM inertia starts to affect the jet advance, two shocks may be 
formed, one moving backwards in the jet, the so-called reverse shock, and another 
one moving forward, the so-called forward or bow shock. Under these conditions, 
Fermi I type acceleration mechanism seems the most reasonable option, although 
high diffusive and convective rates in the downstream regions of the forward/reverse 
shocks could prevent efficient acceleration to occur. It might be the case as well that 
hydrodynamical instabilities distorted the jet and mi xed jet matter with the ISM 
without forming strong shocks (e.g. Heinz & Sunyaev^Sl). 



3.4. Radiative processes 

In the jet base, depending on the dominant conditions, the relevant lep tonic radia- 
tive mechanisms could be synchrotron emission (e.g. Markoff et al.^21), relativistic 
Bremsstrahlung from electrons interacting with jet ions (Bosch- Ramon et al. ^j), 
SSC (e.g. Bosch- Ramon & Paredes^) and IC with corona and/or disk photons (e.g. 
Romero et al. Georganopoulos et al. ^j) . Regarding hadronic processes, there 
are several radiative mechanisms that could produce gamma-rays, neutrinos and, as 
a by-product, low energy emission from secondary pairs. Two of these mechanisms 
are the collisions of relativistic protons with ions (pp) in the jet, and interactions 
between jet relativistic protons and X-ray ph oton s (photo-meson) from the disk, 
the corona or the jet itself (e.g. Levinson et al. Aharonian et al. Romero & 
Vila -who also account for proton synchrotron emission-). These relativistic pro- 
ton collisions with ions and photons would produce neutral pions (tt°) that decay to 
gamma-rays, and charged pions (n^) that decay to muons and neutrinos, the former 
decaying then to electron-positron pairs and neutrinos. Another possible hadronic 
mechanism is photo-disintegration, which requires the presence of UHE heavy nuclei 
and a dense field of target photons of large enough energy. This process produces 
lower mass hadrons and gamma-rays. 
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At binary system scales, possible radiative leptonic p roce sses taking pla ce i n 
microquasars are synchrotron emission (e.g. Yuan et al. ^3 Paredes et al. ^j), 
relativistic Bremsstrahlung (e.g. Bosch- Ramon et al.^j), SSC (e.g . Atoy an & Aha- 
ronian Dermer et al. and external IC (e.g. Paredes et al. Q3E1 Kaufman 
Bernado et al. ^21 Georganopoulos et al. ^2 Dermer et al. EOl Khangulyan et al. 

At these spatial scales, jet proton collisions with target nuclei of the stellar wind 
(e.g. Romero et al. EH, Aharonian et al. E3) seem to be the most efficient hadronic 
process. As noted abov e, this mechanis m wo uld lead as well to neutrino production 
(e.g. Romero et al. Aharonian et al. ^J). Other hadronic processes, which were 
also discussed in the literature, are ph oto- meson production (e.g. Aharonian 
and photo-disintegration (e.g. Bednarek 1^3) . 

The emission at larger scales are commonly characterized by synchrotron radi- 
ation. At higher energies, stellar IC scattering is quite inefficient because the large 
distances to the companion star and the subsequent dilution of the stellar pho- 
ton field. Nevertheless, for powerful ejections, SSC could still be significant (e.g. 
Atoyan & Aharonian Regarding the particle energy distribution, its evolution 
is likely dominated by convective and adiabatic energy losses (van der Laan 12 ^). At 
the termination of the jet, in case particle acceleration and confinement were effi- 
cient, synchrotron, relativistic Bremsstrahlung and IC radiation could be produced 
and even detected for galactic sources (e.g. Aharonian & Atoyan ^\ Bordas et al. 
ISHl). Hadronic acceleration could take place as well, which could lead to gamma- 
ray production (e.g . He inz & Sunyaev and secondary leptonic emission (e.g. 
Bosch- Ramon et al. D3H1) 

Regarding the variability of the emission discussed above, several factors are rel- 
evant: injection could change via variations in the accelerator (e.g. injection power, 
injection spectrum); the target densities could vary via changes in the magnetic, 
photon and matter fields; geometry changes, due to e.g. orbital motion or jet pre- 
cession, could affect anisotropic gamma-gamma absorption and IC scattering or the 
level of radiation Dopplcr boosting. Thus, the timescales of the variability could be 
linked to the injection mechanism, the radiative cooling, particle escape, and the 
orbital motion, depending on which mechanism plays a major role. Depending on 
the emitter location and size, some kinds of variability will not play a role because 
they will be smeared out. We come back to this issue for a more detailed discussion 
in Sect, g] 



3.5. Photon-photon absorption and electromagnetic cascades 

Close to the compact object, either in the inner accretion disk, in the corona- like 
region, or in the jet base, quite extreme conditions should be present. Magnetic 
fields are likely high, and the same applies to the present photon fields, which may 
be dominated by thermal radiation from accretion peaking at UV/X-rays. Dense 
UV/X-ray photon fields imply that ~ GeV gamma-ray absorption close to the 
compact object will be very high. The strong magnetic field may suppress efficient 
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electromagnetic cascading, although the occurrence of cascades (e.g. Akharonian et 
al. ^j) cannot be discarded. At binary system scales, the photon-photon opacities 
of the stellar photon field for gamma-rays are high in massive systems. Although 
high-mass stars can have quite large magnetic fields in their surfaces, up to 10 3 G 
for a young O star (e.g. Donati et al.^j), electromagnetic cascades may still develop 
in the system since the magnetic field strength at several stellar radii is not well 
known. Several authors have studied th e ph oton-photon absorption effe cts in the 
VHE lightcurve (e.g. Protheroe & Stanev^l Moskalenko & Karakula^Sl Bed nare k 
Boettcher & Dermer^ Dubus^l Khangulyan et al. ^ Reynoso et al. H2QJ). 
The effects of absorption on the radiation variability are important not only because 
the photon column density changes along the orbit, but also due to the angular 
dependence of the cross section and the low-energy threshold of the pair creation 
process. The interaction angle between gamma-rays and stellar photons changes 
with the orbital phase. Several studies have been done in the recent years studying 
the impact of cascading in the gamma-ray spectrum of microquasars. Some authors 
have assumed that particles get deflected aft er cr eation, using a three-dimensional 
code to compute cascading, (e.g. Bednarek ^™ ), whereas others have computed 
one-dimension al ca scades in the direction to the observer (e.g. Aharonian et al. 
Orellana et al.^21 Khangulyan et al. ^). In case the magnetic field is high enough, 
synchrotron secondary radiation produced in the sys tem will be significa nt in the 
radio and X-ray domains (e.g. Khangulyan et al. Bosch- Ramon et al. ^03 |) j n 
the next section, we discuss further this issue. 



4. Theoretical interpretation of observations 

To interpret observations, one can use some analogies comparing one specific object 
with the general properties of a population of sources. Unfortunately, neither our 
knowledge on the physics of microquasars emitting at TeV, nor the number of these 
sources, are sufficient to follow this strategy. Another approach, which seems to us 
more suitable to study the complexity of the behavior of these objects at very high 
energies, is to start from basic physical processes to understand the observations in 
a rough but sound way. It implies to check, with the help of data and the available 
knowledge on each particular source, which are the most reasonable physical condi- 
tions, and processes, that lead to the observed VHE emission. To perform such an 
analysis, it is required to constrain the conditions in which particle acceleration is 
possible up to the observed energies, to know which mechanism is behind gamma- 
ray production (either leptonic -IC- or hadronic -pp, P7 and nuclei disintegration-), 
and the impact of photon-photon absorption (and the subsequent energy release 
channel -synchrotron or IC-). In addition, it is important to explore the impact of 
particle transport (advection, diffusion or particle escape) on the final spectrum. In 
the following, each one of these elements is considered individually. It is worth to 
remark that in reality the mentioned processes will be likely coupled, and a very 
complex outcome can be expected, as shown in Sects. 14.2.11 14.3.11 and 14.3.21 when 



July 23, 2009 18:30 WSPC/INSTRUCTION FILE rev 



12 V. Bosch-Ramon & D. Khangulyan 

focusing in the leptonic scenario. 



4.1. Non-thermal processes 

In this section, we explore those processes that are relevant to understand the non- 
thermal emission produced at very high energies for the expected conditions given 
in microquasars. We do not focus particularly in the jet but also in its environment. 
In particular, we consider particle acceleration, particle energy losses and transport, 
and photon-photon absorption and electromagnetic cascading. 



4.1.1. Particle acceleration 

Our aim is not to model accurately what is observed, but to derive basic and solid 
constraints from the available data. The first step is to set a necessary condition for 
acceleration of a charg ed pa rt icle of a certain energy to occur, which is given by the 
Hillas criterium (Hillas u ). This consists on the fact that particles can only be ac- 
celerated if their Larmor radius (tl = E/qB a ; where B a is the accelerator magnetic 
field, and q and E are the charge and energy of the particle, respectively) is smaller 
than the accelerator size (i a ), since otherwise they will escape the accelerator. This 
limits the highest achievable energy to: 

E < qBJ a , (1) 

which can be rewritten in the following form: 

E < 30 ? -B a ,G hi TeV , (2) 
e 

where e is the electron charge, B a ,Q is the accelerator magnetic field (B a ) in Gauss 
units, and h% — ^ a /10 11 cm. Nevertheless, to determine whether particles can be ac- 
celerated up to a certain energy, the specific acceleration and energy loss (or particle 
escape) mechanisms are to be known: t acc = £ coo i/ csc . In general, the acceleration 
timescale can be expressed as: 

i acc = - 0.1 V -E TcV B^ G s , (3) 
c q ' 

where 77 is a dimensionless phenomenological parameter (or function) represent- 
ing the acceleration efficiency, and £/Tev is the particle energy in TeV units. The 
particular case of -q — 1 corresponds to the shortest possible acceleration time in- 
dependently of the acceleration mechanism. Another instance for rj can be given 
for the case of non-relativistic diffusive shock acceleration (plane shoc k with weak 
magnetic field, in the test particle approximation) (e.g. Protheroe ^^ : 

V = 2n-^(—)\ (4) 

where w s h is the shock velocity, and D is the diffusion coefficient (I?Bohm in the 
Bohm limit). For v s h — 3 x 10 9 cm s _1 and D = I?Bohm, V ~ 10 3 . 
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As an example of the importance of using acceleration constraints, we show in 
Fig. [3] a 2-dimensional d a — _B a map for a compact high-mass microquasar with 
stellar luminosity L* = 10 39 erg s" 1 , d, = 3 x 10 12 cm and fcT* w 3 eV. cL is the 



distance from the accelerator to the star. As done by Khangulyan et al. 1^21, it j s 
possible to restrict d a and r\. In the case of LS 5039, for instance, for reasonable 
T) > 10, and the detected 30 TeV photon energies, the accelerator should be outside 
the system, i.e. d a > 2 x 10 12 cm. 
Particle transport 

Without focusing on any particular acceleration mechanism, we note that, either 
the maximum energy is limited by cooling, which implies that particles radiate most 
of their energy before escaping the accelerator, or their escape stops acceleration. 
In the latter case, the emitter itself can be considered larger than the acceleration 
region. The escape time (t esc ) can be characterized, in the accelerator as well as in 
the whole emitter, by the minimum among the diffusion (tdig) and the advection 
(i adv ) times: 

t csc = min(i diff , i adv ) , (5) 
which can be expressed either as: 

for 1-dimensional diffusion in the accelerator/emitter, where I (l±2 — Z/10 12 cm) is 
the (accelerator/emitter) size covered by diffusion, and Bq is the emitter magnetic 




Fig. 3. 2-dimensional d — B map that shows the maximum achievable energy for n = 10 for 
different d a and £? a values. The adopted parameter values are L t = 10 39 erg s — 1 , = 3 X 10 12 cm 
and kT t m 3 eV. 
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field; or as: 



-lO^V^s, (7) 



Vadv 

for advection, where Vadv is the advection speed (V% = Vldv/10 8 cm s _1 ). 

Under the impact of cooling, the typical distance up to which particles can 
propagate is: 

1/2 



, 1n l pl/2 R -l/2 1/2 ( _D_ \ 

l dm = 10 E TcV B G t cooi — 

V -^Bohm / 

under diffusive transport, and: 



cm, (8) 



*adv = 10 10 ( 1Al ^ adv . ) tcooi cm , (9) 



10 10 cm/ S/ 

under advective transport. i coo i is the cooling timescale (in s) of the dominant loss 
mechanism. We discuss in the following different cooling processes that may be 
relevant in the microquasar context. 

4.1.2. Radiative processes 
Synchrotron 

In the presence of a disordered magnetic field, electrons radiate via synchrotron 
emission. The synchrotron cooling timescale is approximately: 

4 



t sy «4x 10 2 (-j (m/m c ) 4 B- 2 G E T ^ s , (10) 

where m and m c are the mass of the particle and the electron, respectively. This 
gives the maximum particle energy: 

/ \ 3/<2 

£ max «60m (m/m c ) 2 (r / B a , G )- 1 / 2 TeV. (11) 

This process is only relevant for leptons unless the magnetic field is very strong, in 
which case hadronic synchrotron may become efficient. 

To see what may happen in a microquasar regarding E max under synchrotron 
cooling, we take B q = 10 6 G as a reasonable value for the jet base magnetic 
field (at ~ 50i?sch ~ 10 s cm from the compact object), obtaining: E max w 
0.06 77- 1/2 TeV. If we assume a distance dependence of B of the form 1/Z and 
locate the accelerator at 10 cm, we obtain: £ max w 6rT 1/2 TeV, i.e. in such a 
situation, it seems unlikely to produce VHE leptons in the jet even at the system 
scales. 

At spatial scales similar or smaller than the binary system, it seems reasonable to 
expect large magnetic fields related either to the jet, to the accretion disk, or to the 
companion star. It implies that synchrotron emission can be an efficient radiation 
process, and electrons may release most of their energy via synchrotron emission. 
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At larger scales, if no significant magnetic field enhancement takes place, efficiencies 
will decrease strongly. 

It is worthy using Eqs. JSJ) and ([9]) to calculate the impact of synchrotron losses 
on the propagation of electrons via diffusion and advection, since they tell how far 
relativistic electrons can go without assuming additional acceleration and neglecting 
other sources of losses. Under the next parameter choice, e.g. B ~ 1 G, T4dv = 
10 10 cm s _1 , D = I?Bohm: advection is the most effective transport mechanism, and 
TeV particles may reach distances of ~ 3 x 10 12 cm from their injection point. 

Inverse Compton 

Under the radiation field of the primary star in a high-mass microquasar, or 
in the case of strong accretion disk/corona emission, IC scattering is to be consid- 
ered and may limit particle acceleration. Synchrotron self-Compton could become a 
dominant process, although it is not treated here, since it would require significantly 
large magnetic fields, making acceleration up to very high energies unlikely due to 
strong synchrotron/ Thomson IC energy losses. In addition, once the electron enters 
in the KN IC regime, the interaction efficiency reduces strongly and synchrotron 
cooling becomes dominant, making SSC not very efficient in the TeV range under 
UV photon fields. We note that, due to the same reason as in the synchrotron case, 
IC losses are only significant for lepton^EI. 

For a Planckian distribution of target photons with temperature T, the IC energy 
loss rate can be approximated, with an accuracy of less than a 3%, by: 



where T mcc = kT/m e c 2 . 

Regarding the particle maximum achievable energy, at the energies in which we 
are interested in this work, IC scattering proceeds in the Klein- Nishina (KN) regime 
(7 3> 1/kT in m c c 2 units). In such simple power-law fit for the cooling 

time can be used for a black-body type of target photon distribution: 



where R and i?* are the distance to the origin and the radius of the source of target 
photons, and T4 = T»/10 4 K. Eq. (|13|) can be expressed in a more convenient form 
in case of a hot and massive primary star as the dominant source of target photons: 



where w — lOOwioo erg cm 3 is the target photon field energy density, and T* the 
stellar temperature. This expression gives a reasonable agreement for the IC cooling 



c For hadrons, other cooling processes, mainly photo-meson production, discussed further in the 
text, will overcome hadronic IC scattering. 




(12) 




(13) 




(14) 
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time, within a factor < 2, in the relevant energy range: 

0A TeV (sISk) 1 < * < 3 x 103 TeV (^k) 1 ' (15) 

The corresponding value of the maximum energy is: 

E max « 4 x 10 10 [B a GV'^wo?- 3 TeV , (16) 

which shows that KN IC limits particle acceleration much less than synchrotron ra- 
diation, and can be hardly dominant at the maximum particle energy for reasonable 
magnetic fields in any region of the microquasar. 

To estimate the impact of IC losses on the propagation of electrons via diffusion 
and advection, Eqs. ([8]) and Q can be used as well. Like in the case of synchrotron 
losses, for the conditions B ~ 1 G, \4dv = 10 10 cm s _1 and D — -Deohm, and IC 
loss dominance, advection is again the most efficient transport mechanism, under 
which TeV particles could reach distances up to ~ 10 12 cm. 

Proton-proton interactions 

As mentioned in the previous section, purely hadronic processes like pp inter- 
actions have been discussed in the past in the context of microquasars. We con- 
sider them here as well, since they may be r elevant in some cases. The energy loss 
timescale for pp collisions (Kelner et al. USUI) j s; 

10 6 

*pp ~ s ' ( 17 ) 

n 9 

where ng — n t /10 9 cm~ 3 is the target density. The energy threshold of this process 
in the reference frame of the interaction center of masses is the pion rest mass, 
ss 140 MeV. From Eqs. (|18p and (|17[) . the maximum particle energy can be derived: 

E max « 10 7 ^ TeV . (18) 

V n 9 

Given the long cooling timescale, the maximum energy will be likely limited by the 
accelerator size. 

Defining L p as the luminosity injected in the form of relativistic protons, the 
luminosities in gamma-rays, neutrinos and seconda ry e lectron-positron pairs are, 
with differences of about a factor of 2 (Kelner et al. ^ ^he next ones: 

L 1 k min(l,te SC /*pp)cppZp, (19) 

where c pp is the energy transfer efficiency from relativistic protons to secondary 
particles (~ 10%). In the context of high-mass microquasars, a reasonable lower 
limit for t esc is the wind advection time, i.e. the time required for the stellar wind 
to cross the orbital radius (R rb/V w ; where i? or b ~ 10 12 — 10 13 cm and V w ~ 
1 — 3 x 10 8 cm s _1 ), t csc < 10 4 s; and a very lower limit is set by the speed of light, 
i.e. tesc > 10 2 s. All this, plus adopting a M = 10 -6 M yr _1 (typical for O stars), 
yields: 

L 7 ~ 10" 5 - 10~ 3 L P . 
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Larger efficiencies cannot be discarded in some specific cases, like density enhance- 
ments in the wind or even in the jet itself via, e.g., shocks. Once relativistic protons 
leave the binary system, the expected density decreases in the stellar wind as well 
as in the jet, making pp collisions negligible. 
Photo-meson production 

Among different hadronic processes, photo-meson production (Kelner et alJ ^^ j) 
is worthy also to be considered. The energy threshold for this process is: 

E th P7 = m p c 2 e th P7 /2e = (5 x I0 4 TeV) (T^ 1 , (20) 

where m p is the proton mass and e and eth p 7 ~ 140 MeV are the energy of the 
target photon in the laboratory and the hadron rest frames, respectively. The loss 
rate is given by: 

10 18 

N x 

where: 



*P7~— S > ( 21 ) 



A^x « -^r c « 2 x 10 14 L 3S T^R^ 2 cm" 3 . (22) 

1/38 is the star luminosity L*/10 38 erg s _1 , R12 is the distance to the star, and 
i?/10 12 cm. Eq. (f2"Tj) can be rewritten as: 

Vy - 10 4 i 38 1 T 4 i? 2 2 s . (23) 

With this cooling timescale, the corresponding maximum energy is: 

£ max ~ IOV^g^TA TeV . (24) 

As for pp interactions, the long cooling timescales of photo-meson production 
imply that the maximum energy is in fact limited by the accelerator size. We note 
that only in case a substantial part of the energy in accelerated protons were > 
Eth pm ~ 10 4 TeV, under an UV photon field, and the fastest protons escaped at 
the speed of light, photo-meson production could reach efficiencies of L 1 ~ 10~ 3 LpE| 
for E p > -E>th p7- The energy transfer efficiency for this process is c P7 = 0.1. In 
the inner regions of the jet, close to the accretion disk and corona photon fields, the 
radiation energy density may be high enough to get even higher efficiencies. The 
larger target photon energy would imply a reduced threshold energy. Nevertheless, 
the constraint imposed by Eq. @ is quite restrictive and could prevent hadrons 
from reaching energies > -Eth p 7 for reasonable B a and Z a values. Namely, comparing 
Eqs. (|24p. (|20p and two conditions are obtained: 

B G T A l n > 2 x 10 3 , (25) 

for basic physical conditions; and 

V < B aiG TiR 2 12 /2L 3S . (26) 



d At the involved proton energies, the stellar wind can hardly confine the particles, thereby we 
adopt the speed of light as escape velocity. 
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for the acceleration efficiency. 

As in the case of pp collisions, the strong dilution of the photon field far from 
its source prevents photo-meson production from being efficient outside the binary 
system. 

Photo-disintegration 

If ultra relativistic heavy nuclei are present, they can suffer photon disintegration 
under the ambient photon field. The expression for the threshold energy is similar 
to that presented in Eq. (|20jl : 

E t h P d = m N c 2 e t h P d/2e , (27) 

where ton is the mass of the nucleus and e t h pd = 8 MeV. Effectively, since ton can 
be up to ~ 100 to p , Eth pd could be > Eth pj- 

The disintegration of the nuclei has as typical timescale: 

t pd ~ 10 3 L^T A R 2 12 s . (28) 

For simplicity, we have assumed that the mass of the nucleus is its charge times the 
mass of the proton, which slightly overestimates the efficiency within a factor of ~ 2. 
In addition, a slow dependence on the mass of the nuclei in the cross section has 
been neglected. Photo-disintegration would stop the acceleration of heavy nuclei at 
energies: 

Bmax ~ W^B^L^T^ TeV . (29) 

As in the case of pp collisions and photo-meson production, the accelerator size in- 
stead of photo-disintegration cooling will probably stop acceleration of heavy nuclei. 
Taking into account that the energy transfer efficiency is c P7 ~ 0.01, and adopting 
the speed of light to derive a lower-limit for t csc , photo-disintegration can yield in 
gamma-rays L 1 ~ I0 _3 Ln, where Ln is the energy stored in heavy nuclei with 
E > E^ pd- Hillas constraint allows energies ton/to p times larger for nuclei than 
for protons, although as noted above E t h P d can be > Eth p-JE Like in the case of 
photo-meson production, it is unclear whether enough energy can be in the form 
of ultra relativistic heavy nuclei above the threshold to produce a significant signal 
via photo-disintegration. 

Photon-photon absorption and secondary radiation 
The presence of the star, a powerful source of UV photons, is very relevant 
to understand TeV emission from compact sources. A massive and hot star is an 
excellent provider of target photons for photon-photon absorption. The anisotropy 
of the stellar photon field and the relative position of the system with respect to 
the observer makes the exploration of the observational impact of photon-photon 
absorption a non-trivial subject. In Fig. [4j we show a 2-dimensional map of the 
opacity coefficient (r) for different energies of the incoming photon (Ey) emitted 



c Actually, similar constraints to those derived for photo-meson production regarding the interac- 
tion physical conditions and the acceleration efficiency apply here as well. 
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near the compact object, and different directions starting with the one pointing 
away from the star. Here, 9 is the angle with the line joining the TeV emitter and 
the star. The point-like approximation for the target photon field has been used, 
which implies that r oc 1/d. The used parameter values are L* = 10 39 erg s _1 , d, = 
3 x 10 12 cm and fcX 1 , = 3 eV. It is worthy noting that the point-like approximation 
for the target photon emitter works in a wide range of situations, but there are cases 
relevant for compact TeV sources in which is necessary to account for the finite size 
of the star (see also Dubus ").In Fig.0 we show a 2-dimensional map of the ratio 
Tp/Vf (i.e. point-like versus finite size r) for different distances to the star (d) and 
values of 9. The radius of the star has been taken i?» = 7 x 10 11 cm. The energy of 
the incoming photon has been taken 1 TeV. The deep blue/black regions in this map 
represent the cases when the point-like approximation gives almost no absorption, 
i.e. when the emitter points roughly away from the star, or when the photons in the 
finite size approximation collide with the star surface giving infinite opacity. 

The importance of cascading has been already stated above. It is worth now to 
study the possibility of cascade development in the surroundings of massive and hot 
stars. There are two extreme situations. Either the magnetic field energy density is 
much lower compared with the radiation energy density, and in such a case pure 
cascade will develop, or the magnetic field energy density is a significant fraction, or 
above, the radiation one, and then the electron energy will be mostly released via 
synchrotron radiation and single-scattering IC. If the magnetic field is low enough 
(and therefore suitable for cascading to occur), very high-energy particles will be 
only slightly deflected before radiating their energy via IC. This allows us for a 
1-dimensional approach to compute cascading. For the same parameter values as 
those presented above, we have performed a 1-dimensional electromagnetic cascade 
simulation, the result of which is shown in Fig. [5] The injected particle spectrum 
was a power-law of photon index 2. 

To give an idea of the magnetic field importance, next formula shows the critical 
value of the magnetic field that allows cascading to occur for TeV electrons: 



In fact, B c is valid for the 1-dimensional cascade case. If electrons suffer strong 
deflection in the ambient magnetic field, the longer time required by particles to 
escape the region of dense photon field will increase the synchrotron outcome with 
respect to the IC one. This happens since photons convert to electrons more times 
before escaping, releasing at the end more energy in form of synchrotron and low 
energy IC radiation. From all this, and the fact that magnetic fields of hundreds of 
G are typical in the surface of OB stars, we conclude that effective electromagnetic 
cascading is unlikely in the environment of high-mass systems. For the same condi- 
tions as those taken for Fig. 6, plus a stellar surface magnetic field of 100 G and a 
primary gamma-ray injection luminosity of 3 x 10 35 erg s _1 , we show the secondary 
pair spectral energy distribution (SED) in Fig. [7J We note the moderately high 




(30) 



July 23, 2009 18:30 WSPC/INSTRUCTION FILE rev 



20 V. Bosch-Ramon & D. Khangulyan 



radio and X-ray luminosities (see also Bosch-Ramon et al. 



4.2. Modeling non-thermal emission in microquasars: the case of 
Cygnus X-l 

In Sect. 14.11 we have reviewed the main radiative processes that could take place 
under typical conditions in microquasars. After roughly estimating the efficiency 
of different leptonic and hadronic processes, it is our aim now to focus on the 
mechanism that, to our understanding, is the most likely to produce the VHE 
radiation observed in some of these sources. From the point of view of efficiency, 
IC scattering is clearly a good candidate in this sense. In addition, the fact that a 
hot and massive star is present in all these systems detected at TeV energies makes 
IC scattering very attractiv^. Finally, the pattern of variability of the spectrum 
and flux along the orbit found in LS 5039 can be, together with photon-photon 



It does not discard low-mass microquasars as VHE emitters. Actually, from the observational 
point of view, hints of TeV emission from GRS 191 5—10 5 were found using imaging C heren kov 
techniques by HEGRA (Aharonian &: Heinzelmann EBB ) and Whipple (Rovero et al. ESI). I n 
addition, a n umb er of theoretical wor ks h ave proposed t hese sources as VHE emitters (e.g. Atoyan 
& Aharonian E3 Bosch- Ramon et al.ES Romero et al.ESl). 

100 
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Fig. 4. 2-dimensional 9 — B gamma map of the opacity coefficient. The adopted parameter values 
are L* = 10 39 erg s _1 , d,=3x 10 12 cm and kT, = 3 eV. 
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absorption, explained due to the angular dependence of the cross section of IC 
scattering an d the anisotropic nature of the stellar photon field (e.g. Khangulyan 



To carry out a simple but thorough treatment of the processes relevant to VHE 
emission in microquasars, we have simplified the scenario as much as possible. We 
show a schema of this scenario in Fig. [8] We consider the jet as a weakly relativistic 
convective flow without going into the details of its physical nature, i.e. whether 
it is a magneto- or a hydrodynamical plasma, or even a pointing flux dominated 
flow. This outflow is ejected perpendicular to the orbital plane bearing a disordered 
magnetic field attached to it. Efficient acceleration of non-thermal particles can 
occur at a certain location in the jet (Zq from the compact object, do from the 
star); such acceleration regions are treated as point-like, and their position in the 
jet, Z , could change. For simplicity, we adopt a model in which acceleration takes 
place only in one point. 

Because of the strong uncertainties affecting the physics of microquasar jets, we 
do not treat but very generally the acceleration processes that may take place in 
there. The injection electron spectrum adopted here is phenomenological: a power 
law of exponent -2 plus an exponential high-energy cutoff. The accelerator sets up 
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Fig. 5. 2-dimensional 9 — d map of the ratio r p _i/rfj n for an incoming photon of 1 TeV. The 
adopted parameters are the same as those of Fig. [4] plus the stellar radius R* = 7 X 10 11 cm. The 
shaded area to the left of the plot corresponds to a region of opacities ~ 0. 
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the initial conditions of the emitter, which is treated as a 1-dimensional region 
where particles are transported by jet convection. Adiabatic losses have been ne- 
glected here. The magnetic field depends on Z like B = B a (Zo/Z). To model the 




Photon energy (GeV) 



Fig. 6. 1-dimensional electromagnetic cascade. Note the difference of about a factor of 10 between 
the 1 TeV and the 10 GeV fluxes. The adopted parameters are those of Fig. [4] 
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Fig. 7. Computed SED of the synchrotron and IC emission from the pairs created in a compact 
binary. The TeV emitter location is at d = 3 X 10 12 cm. The adopted parameters are those of 
Fig. [4] plus a stellar surface magnetic field of 100 G and a primary gamma-ray injection luminosity 
of 3 X 10 35 erg s" 1 . 
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non-thermal processes included in our scenario, it is necessary to consider the pres- 
ence of different ingredients: the magnetic field, possibly present radiation fields, 
and the material in which the emitting particles may be embedded. As discussed 
above, we will deal here with synchrotron and stellar IC emission, since we restrict 
ourselves to high mass systems and to the most relevant processes to produce or 
affect VHE radiation. Nevertheless, other mechanisms cannot be discarded. We take 
into account the impact of photon photon absorption on the VHE spectrum. We re- 
call the importance of the geometry of IC scattering and photon photon absorption 
for the spectra and lightcurves when the star is the main source of target ph oton s. 
For further mathematical details of the model, we refer to Khangulyan et al. ^1 

In the following, we apply this model to the microquasar Cygnus X-l 
(Sect. 14.2. and to the other two X-ray binaries with extended radio emission, 
LS 5039 and LS I +61 303 (Sect. 0}. 

4.2.1. Application to Cygnus X-l 

Cygnus X-l is a HMXB of uncontroverted accreting nature with relativistic radio 
jets (Stirling et al. ^^ j). The compact object is a black- hole of ~ 20 M , the primary 
is an 09.7Iab star of ~ 40 M©, and the system is located at 2.1 kpc (Ziolkowski 
EH). The orbit is thought to be circular, with an inclinatio n i ~ 30°, orbital radius 
of « 3 x 10 cm and period of 5.6 days (Gies & Bolton EH) xhc 

primary star 

presents a luminosity of « 10 39 erg s _1 and a temperature of « 30000 K. An 





Fig. 8. Schema of the system. 
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evidence of detection from this source at the 4.1 a level has been found by MAGIC 
in the TeV range during a tr ansi ent event that may be correlated with the hard 
X-ray behavior (Albert et al. US). At present, the origin of the VHE emission is 
uncle ar. In the context of the hadronic scenario, Romero et al. I^H and Orellana 
et al. computed the emission for a microquasar with the characteristics of 
Cygnus X-l. Regardless the hadronic or leptonic origin of the gamma-ray radiation, 
electromagnetic cascading and / or secondary s ynch rotron emission should occur in 
the system (e.g. BednarekEH, Orellana et al. ISH, Bosch-Ramon et al. I ^Q^ I). 
On the spectral energy distribution and orbital variability 
It is worthy to see the aspect of the multiwavelength non-thermal emission in- 
ferred from the VHE radiation detected from the source. This is shown in Fig. [9] 
We have located the accelerator/emitter at Zq = 3 x 10 12 cm. The magnetic field is 
~ G. The synchrotron X-ray flux is well bel ow t he observed level, of likely thermal 
(comptonized) origin (Sunyaev & TruemperE^). This SED has been calculated for 
an emitter located in the borders of the system, in which opacities are moderate. 
However, deeper, closer to the compact object, the attenuation factor grows orders 
of magnitude. 

For Cygnus X-l, the spectrum of the emission from ~ 0.1 — 1 TeV electrons is 
available for a narrow orbital phase range around 4> = 0-9, i.e. this is not averaged 
for wide ranges of the orbit. In addition, these electrons have the shortest lifetimes 
under IC cooling, implying that they radiate most of their energy before propagating 
significantly. These two facts mean that the observed radiation was produced under 
similar conditions. In addition, the stellar photon field is very dense in the compact 
object surroundings, implying very large opacities in almost all the directions. All 
this, plus the known orbital observer /system geometry at the observation phase, 
allows us, and makes interesting, to estimate the absorbed energy via photon-photon 
absorption depending on the distance between the emitter and the compact object. 
The calculation of the energy processed in this system via photon-photon absorption 
can give a model-independent constraint on the the location of the TeV emitter. 

In Fig. [lOl we show the luminosity divided by 47rdg (cIq : distance to the Earth) 
of the secondary pairs injected in the system as a function of the distance to the 
compact object. This is computed calculating first the deabsorbed VHE spectrum 
and flux given the distance to the star and the observed spectrum and flux (Albert 
et al. 15 ). In our calculations we have not taken into account the effect of cascad- 
ing, although the likely moderate-to-high magnetic fields present in massive binary 
systems would imply a dominant synchrotron channel suppressing cascading. Since 
in our context a significant fraction of this energy rate is released via synchrotron 
emission, and the X-ra y /sof t gamma-ray luminosities are typically of ~ 10 37 erg s _1 
(e.g. McConnell et al. ^^ j), the emitter location ca n be set, conservatively, to at 
least few times 10 12 cm (see also Bosch- Ramon et al. ^^ ). This is a strong indica- 
tion that the TeV emitter is located in a jet far away from the compact object. It 
is relevant also to note that the secondary synchrotron r adiat ion may also explain 
the observed soft gamma-ray fluxes (e.g. McConnell et al. ^^ j). Another interesting 
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result is the absorbed luminosity at distances ~ 10 13 cm, indicating that there can 
be significant injection of relativistic pairs at distances where their r adio emission 
may be resolved using VLBI techniques (see also Bosch- Ramon et al. 

Regarding the time evolution of the observed radiation, the flaring nature of 
the VHE emission points to intrinsic changes of the emitter properties, instead of 
geometric effects or target density variations, as the origin of the variability. 



4.3. TeV binaries with extended radio emission: LS 5039 and 
LS I +61 303 

4.3.1. LS 5039 

LS 5039 is an HMXB (Motch et al.EHI) loc ated a t 2.5±0.5 kpc (Casares et al.Qlj). 
The source presents radio jets (Paredes et alJ ^^ j), s how s X-ray variability possibly 
associated to the orbital motion (Bosch- Ramon et al. ~~h and has been detected at 
very high-energy (VHE) gamma-rays (Aharonian et al.L^), w hich virtually confirms 
its association with a CGRO/EGRET source (Paredes et al. ^*). Interestingly, the 

n~7i\ 

TeV emission varies with the orbital period (Aharonian et al.L^. The most recent 
orbital ephemeris of the system were obtained by Casares et al. The orbital 
period is 3.9060 days, the eccentricity is moderate e = 0.35±0.04, and the inclination 
angle is not well constrained, i » 15°-60°. The orbital semi-major axis of the system 
is w 2 x 10 12 cm. The nature of the compact object is still uncertain. Casares et. 




E, TeV 



Fig. 9. Computed synchrotron and IC SED for Cygnus X-l, together with the observed VHE 
SED, is shown. The accelerator/emitter location has been set to = 3 X 10 12 cm. The total 
contribution from the jet and the counter-jet is shown. 
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al. ^Hl suggested that it may be a black holeQ, although there is an on-going debate 
on this issue, and some authors have proposed that LS 5039 is in fact a young 
non-accreting pulsar (see, e.g., Martocchia et al. see a j so D u bu S 11321) . In this 
regard, the strongest argument would the lack of accretion features in the X-ray 
spectrum of LS 5039. 

The radio emission in LS 5039 is of non-thermal origin, slightly variable at 
month-year timescales (~ 20-3 0%) and extended - ~1 - 1 00 milliarcsec onds (mas)- 
(Martf et al. GM Ribo et al. ESQ Paredes et al. UMH^ mho et al. 122 for an 
exhaustive study of this source, see Ribo^^l), the ~ 60-80% being produced within 
a core of ~ mas. The radio emission, when observ ed a t VLA angular resolution, 
appears unresolved, and optically thin (Marti et al. Q2H) 

In the X-rays, the source shows flux variations by a factor of ~ 2 peaking 
smoothly around phase 0.8 and more sharply at other phases (Bosch-Ramon et al. 



s Assuming pseudo-synchronization of the orbit. 




Fig. 10. Luminosity divided by 4-7rdg of the secondary pairs injected in the system with the 
emitter at different distances from the compact object. This is computed calculating first the 
deabsorbed VHE spectrum and flux given the distance to the star and the observed spectrum and 
flux. The case in which the radiation process is isotropic, like proton-proton collisions may be (solid 
line), and the case in which the radiation process is IC (dot-dashed line) -strongly anisotropic-, 
are shown. The production curves for both cases are also shown. 
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. These peaks were apparently accompanied by spectral hardening. At phases 
~ 0.8, higher and harder emission than in the rest of the orbit has been also observed 
at TeV energies (Aharonian et al. ^), and simultaneous Cha ndra observations have 
ap paren tly shown a similar behavior in X-rays (Horns et al. Q~251 Bosch- Ramon et 
al.™). 

At VHE, the flux and the photon index change periodically, by a factor of 
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Fig. 11. At the top, the lightcurve of the photon flux above 1 TeV of LS 5039 presented by 
Aharonian et al. (2006). Two full phase periods are shown for clarity. The blue solid arrows depict 
periastron and apastron. The thin red dashed lines represent the superior and inferior conjunctions 
of the compact object, and the thick red dashed line depicts the Lomb-Scargle Sine coefficients for 
the period giving the highest Lomb-Scargle power (see Aharonian et al. 2006 for details). At the 
middle, the fitted pure power-law photon index (for energies 0.2 to 5 TeV) versus phase interval of 
width <j> = 0.1 is presented. At the bottom, the power-law normalization (at 1 TeV) versus, phase 
interval of width A0 = 0.1 is shown. 
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~ 4 the former and between ~ 1.9- 3.1 the latter, with the spectrum hardening 
when flux increases (Aharonian et al. ^j). The observed lightcurve and VHE SED 
are shown in Figs. [IT] and 1121 respectively. The maximum of the emission takes 
place around phase 0.8, similarly as it seems to happen at X-rays. Also, sudden 
increases/hardening in the flux and spectru m on hour timescales could have been 
detected at phase ^0.85 (de Naurois et al. ^3) j being similar to what has been 
mentioned above concerning sharp peaks in X-rays. The TeV emission from LS 5039 
has been studied by several authors in th e microquasar (e.g. Paredes e t al. and 
the pulsar context (e.g. Dubus et al.Q2Sl Sierpowska-Bartosik & Torres ESEIUl ), A 
more general approach to study t he VHE emission from LS 5039 has been recently 
carried out by Khangulyan et al. ^1 

In the following, we present the multiwavelength SED and the orbital VHE 
lightcurve computed with the model sketched above for LS 5039. The values given 
to the free parameters of the model are chosen with the intention to show that a 
simple leptonic model can roughly explain the observed features of the emission 
taking into account the TeV data and the phenomenology of the source. At this 
stage, no electromagnetic cascading should be introduced. 

On the spectral energy distribution and orbital variability 
In Fig. 1131 for illustrative purposes, we show the synchrotron and IC SED for 
LS 5039 averaged over the orbital phase interval 0.9 > cj) > 0.45, corresponding to 
the inferior conjunction of the compact object (<fi = 0.72). We use this format for a 
better comparison with the spectral results from observations (Aharonian et al.03), 
shown also in the figure. For the very high energy spectrum in superior conjunction, 
see Khangulyan et al. ^31, fig. 23; at that phase interval, the synchrotron emission 




E(eV) 



Fig. 12. Spectra of LS 5039 for the two broad orbital phase intervals 0.45 < <f> < 0.9 (red circles), 
and rj> < 0.45 and <j> > 0.9 (blue triangles), from Aharonian et al. (2006). The shaded regions 
represent the 1 a confidence bands on the fitted functions (see Aharonian et al. 2006 for details). 
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will be similar to that at inferior conjunction. In Fig. [Ml the VHE lightcurve along 
the orbit is shown. The adopted parameter values in both figures are the following: 
y adv = 0.1 c, Z Q = 10 12 cm, i = 25°, r\ = 10 and B = 0.05 G. The exact values of 
B and V a dv are actually no crucial. In the case of B, its value just needs to fulfill 
acceleration constraints (i.e. to allow electrons to reach high enough energies to 
explain the TeV data), and be low enough such that VHE electrons will lose most 
of their energy via IC in the KN regime. Regarding V at i v , moderately different values 
will render similar results. To reach the observed emission levels, the total injected 
luminosity in rclativistic electrons is L c ss 1.5 x 10 35 erg s _1 . We note that we have 
computed the contribution from both the jet and the counter-jet and summed them 
up. Photon-photon absorption is taken into account. Remarkably, the synchrotron 
part of the spectrum is well below the observed fluxes, indicating that the dominant 
X-rays in this source likely come from a different emitting region. This region should 
be nevertheless physically connected to the TeV emitter, given the similar behavior 
in time of both the X-ray and the TeV radiation. 

Regarding variability, we show in Fig. [14] the orbital VHE lightcurve obtained 
for the same parameter choice adopted to compute the SED. Although the matching 
with observations is not perfect, it is necessary to state that the lightcurve provided 
in Aharonian et al. ^3 j s constructed in bins with relatively small statistical sig- 
nificance, the spectra are simplified to pure power laws (and they could be more 
complicated than just a power law, as shown in Fig. I13p. and includes data points 
from observations of very different epochs. In any case, just playing with propaga- 
tion effects, IC scattering and photon photon absorption, one can already obtain 
very different lightcurves. 

For an extensive discussion on the role of th e d ifferent parameters present in 
the used model, we refer to Khangulyan et al. ^1 Here our purpose is just to 
note that the leptonic scenario, accounting for changes in the interaction geometry 
due to the orbital motion of the system, can reproduce quite nicely the observed 
spectra for different phases. It is worthy mentioning that, to reach the observed 
maximum photon energies, Zq is to be > 10 12 cm unless the acceleration efficiency 
is close to the maximum one, i.e. r\ < 1(0. In addition, the magnetic field in the 
accelerator /emitter is restricted to a relatively narrow band around 0.1 G. 

4.3.2. LS I +61 303 

LS I +61 303 is a quite eccentric (e = 0.72) HMXB, located at ~ 2 kpc, with an or- 
bital semi-major axis of ~ 6x 10 12 cm, a Be primary star of luminosity w 10 38 erg s _1 
and temp eratu re w 28000 K, and an orbi tal p eriod of 26.4960 days (Hutchings & 
CramptonQSl Gregory 11221 Casares et alJ ^^ j). The inclina tion of th e system is not 
well constraint, being in the range i = 15 — 60°. Massi et al. (| 134 | 135 [j d e t ec t ec i radio 



h In the case of a deep emitter, the m agnetic field in the system should be well below several G, 
as noted by Bosch-Ramon et al. EZ] 
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jets, with apparently relativistic motion, at ~ 100 milliarcsecond scales and classi- 
fied LS I +61 303 as a microquasar. Otherw ise, t he source does not s how signatures 
of accretion in the X-rays (e.g. Sidoli et al. USE Paredes et al. Q2Zl an d references 
therein), which are of likely non-thermal origin. This fact, altogether with other 
argu ments -mainly related to the extended radio emission morphology (see Dhawan 
et al. have let several authors to p ut fo rward a non-accreting pulsar as the 

co mpac t object in the system (e.g. DubusH2Hl Chernyakova et al. 11301 Dhaw an et 
gj 1 138 |)^ w hi c h had been proposed for the first time by Maraschi & Treves Nev- 
ertheless, the microquasar scenario cannot be discarded, and it may indeed explain 
some of the i nferre d properties of LS I +61 303 better than the pulsar scenario (e.g. 
Romero et alJ ^^ j). In fact, the apparently slow and collimated radio structures de- 
tected in this source may be rather difficult to explain in the colliding wind context, 
in which ve ry fa st motions of the shocked pulsar wind could be expected (see Bogo- 
valov et al. ^2T). We recall that LS I +61 303 has been detected in the TeV range 
by MAGIC (Albert et al. with the maximum around the ph ase .6, being not 
detected during periastron passage, at phase 0.23 (Ca sares et al. H22l). This source 
had been also detected by EGRET (Kniffen et al.Q33l). It is worthy noting that the 
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Fig. 13. Computed synchrotron and IC SED of LS 5039 averaged over the orbital phase interval 
0.9 > 4> > 0.45. The HESS data points are also shown (Aharonian et al. 2006). The adopted 
parameter values are V a d v = 0.1c. Zq = 10 12 cm, i = 25°, rj = 10, B = 0.05 G, and L e = 
1.5 X 10 35 erg s _1 . The total contribution from the jet and the counter-jet is shown. 
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lightcurves in the radio, X-rays, and high-energy and v ery h igh-energy gamma-rays 
look somewhat similar (see fig. 3 in Chernyakova et al. ESI and references therein). 
Several authors have adopted different frameworks and mechanisms to explain the 
high energ y radiation from this source (in the microquasar framew ork: e .g. Romero 
et al. 1 1^ 1 -hadronic-, Bosch- Ramon et al. 1 1 ^6 1 -lept onic- , Bednarek -cascading- 
; i n the pulsar framework: e.g. Leahy Dubus ^221 -leptonic-, Cherny akoy a et 
al. 0^01 -hadronic -; re garding neutrino detectability, see Christiansen et al. anc l 
Torres & Halzen EHI) . 

On the spectral energy distribution and orbital variability 
In Fig. HH we show the computed synchrotron and IC SED for LS I +61 303 
at <j> = 0.6. The adopted parameter values are the following: 14dv = 0.1c, Zq = 
10 12 cm, i = 30°, and r\ = 10. In this case, the magnetic field energy density could 
be larger than the radiation one compared to LS 5039, since the VHE spectrum is 
less hard, leaving more room for synchrotron losses (which yield a softer spectrum) . 
The value of B obtained here to reproduce the observed radiation is 0.5 G. The 
total adopted L c is 1.5 x 10 35 erg s _1 . In LS I +61 303, the X-ray levels would 
not be far from the observed ones. Still, for the adopted magnetic field, adjusted 
to explain the TeV spectrum, the fluxes are still too low. Unlike LS 5039, neither 
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Fig. 14. Computed IC orbital lightcurvc of the differential photon flux at 1 TeV for LS 5039. 
The different components, jet, counter-jet and summation of both, are labelled. The adopted 
parameters are the same as those adopted in Fig. 1131 
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the observer/binary system geometry nor the impact of photon-photon absorption 
are enough when understanding the TeV variability, given the distance to, and the 
luminosity of, the primary star. At the phase at which this emission is detected, the 
gamma-ray attenuation is low, and the IC scattering angle changes relatively little. 
In Figs. [16j [TT1 and [TBI the differential photon flux lightcurves at 0.1, 1 and 10 TeV, 
respectively, are shown. In the observed lightcurve (Albert et al. H^j) , the m aximum 
of the emission is at p hases (using the ephemeris of Casares et al. Q2S1 see a j so 
Grundstrom et al. H5QI) different from those predicted by our IC/photon-photon 
absorption model. This fact, plus the similarity of the VHE lightcurve with those at 
other wavelengths, v suggest intrinsic changes of the emitting region. Interestingly, 
the 10 TeV emission in our scenario peaks roughly around periastron passage, due 
to the characteristics of the pair creation and IC cross sections. 

5. Further comments 

We have discussed in this work the different processes that may take place in micro- 
quasars, and explored the possibility that stellar IC scattering is behind the VHE 
radiation from three sources recently detected in the TeV range. We have shown 
also the importance of photon-photon absorption, whose occurrence can be used 
to constrain the emitter properties. In this section, we critically discuss three hot 
topics that are, in our view, very relevant for the kind of source discussed here. The 
first point is related to the possibility, already mentioned and commented above, 
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Fig. 15. The Computed Synchrotron and IC SED for LS I +61 303 at phase (j> = 0.6, together with 
the observed VHE SED, is shown. The adopted parameter values are V a( iv = 0.1 c, Zq = 10 12 cm, 
i = 30°, jj = 10, B = 0.5 G and L c = 1.5 X 10 35 erg s _1 . The total contribution from the jet and 
the counter-jet is shown. 
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that a non-accreting pulsar may power the non-thermal broadband emission in 
LS I +61 303 and LS 5039. The second point, affecting Cygnus X-l as well, is the 
role played by the companion star concerning high-energy processes and the dy- 
namics of the jet. Finally, some comments are made regarding the prospects of the 
study of microquasars with the new/future VHE intrumentation. 



5.1. Evidences for a pulsar nature 

Although we do not consider here that LS I +61 303 and LS 5039 may harbor young 
non-accreting pulsars instead of an accreting black-hole, this possibility cannot be 
discarded. At present, the strongest argument in favor of this is the lack of accretion 
features in the X-ray spectrum. An extension of this argument is the fact that 
the timing properties of the X-ray emission in these two sources, and even the 
radio/X-ray behavior, do not correspond to what is believed a microquasar should 
be. Needless to say, this argument is phenomenological, and relies on a supposed 
regularity of the microquasar behavior at X-rays. The same could be said regarding 
the radio behavior and the radio/X-ray connection. 

Indeed, there is a whole set of phenomenological studies, as already mentioned 
in Sect. [2 majorly grounded on observations but also to some extent on theoretical 
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Fig. 16. The light-curve (differential photon flux) at 100 GeV for LS I +61 303. The jet and 
the counter-jet components are labelled. The adopted parameter values are the same as those in 

Fig. rug 
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models, which try to give an unified microquasar picture. LS 5039 and LS I +61 303 
do not fit in such a picture regarding the points mentioned above. Nevertheless, 
most of the high-mass microquasars, and to a lesser extent several low-mass micro- 
quasars, present some level of discrepancy from the standard picture, and not all the 
sources may have occurring in them the same mechanisms as those comprehended 
in the standard scenario. All this shows that such a picture or framework is a useful 
working hypothesis, but it cannot still be a discrimination tool when trying to find 
out whether or not a source pertains to some class of objects. In the same line, 
sometimes the morphology of the radio emission from these systems is used as an 
argument, again phenomenological, against their microquasar nature. In this regard, 
it is claimed that the extended radio emission does not fit the standard picture of 
microquasar jets. Nevertheless, it is worthy pointing out here that galactic compact 
jets are sometimes detected at the resolution limits of radio VLBI interferometers. 
It makes therefore difficult to argue, based on solid observational grounds, about 
how a microquasar jet should look like. From the theoretical point of view, to define 
how a canonical microquasar jet should be is presently not possible, since we lack 
a complete theory for jet formation and collimation, and there might be more than 
one mechanism operating. 

There are otherwise two good indicators of the presence of a pulsar in the sys- 
tem, i.e. the detection of radio and X-ray pulsations, and the lack of strong X-ray 
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Fig. 17. The same as in Fig. [TBI but at 1 TeV. 
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radiation (as is the case) due to any form of accretion onto a possible neutron star 
surface (in case the compact object were known to be a neutron star). To date, none 
of both questions has been answered. Neither pulsations have been detected, nor 
the problem of the compact object nature has been solved. The lack of pulsations 
may be explained by dense stellar wind ionized material smearing out the pulsed 
signal via free-free absorption, although the powerful pulsar wind would allow the 
observer to see the pulsar without the interference of the stellar wind around inferior 
conjunction. The radiation beam may also point away from us preventing us from 
seeing it. Regarding accretion X-ray bursts, since the masses of the components in 
LS 5039 and LS I +61 303 are not yet properly constraint, we cannot tell whether 
the compact object is a black hole or a neutron star. 

Therefore, the question whether there is an accreting or a non-accreting compact 
object in LS 5039 and LS I +61 303 remains open. Fortunately, because of this 
fact, their complex multiwavelength spectral and temporal behavior, and their TeV 
detections, these two sources are extensively studied nowadays. This will bring for 
sure fruitful and surprising results in the near future. 
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5.2. The relevance of the primary star 

LS 5039, LS I +61 303, and Cygnus X-l contain a massive and hot OB type star 
which embeds the jet/accelerator/emitter with dense matter and photon fields. 
Unlike low-mass microquasars, where the accretion/ejection phenomena could be 
naturally reduced to the disk/jet system, in high-mass microquasars (SS 433 and 
Cygnus X-3 are two additional instances of high-mass microquasars, both being cer- 
tainly peculiar), the strong photon field should play an important role in the radia- 
tion processes via, e.g., stellar IC, photo-meson production and photo-disintegration, 
photon-photon absorption, secondary pair radiation and electromagnetic cascading. 
In addition, the strong stellar wind is likely to have an impact in the radiation pro- 
cesses via dynamical interactions with the jet/accelerator/emitter, providing targets 
for pp radiation, confining relativistic particles, either secondary pairs or protons, 
absorbing part of the radio and X-ray photons produced in the system, and deter- 
mining the medium at large scales (and thereby influencing the processes that take 
place in the termination regions of the microquasar jets). As shown along this re- 
view, numerous studies have been or are being carried out to clarify the importance 
of the primary star in the mentioned processes. 

Beside the fact that the physics to extract from theoretical studies plus present 
and future observations can teach us a lot about jets, acceleration and radiation pro- 
cesses, and stellar winds, it is also important to remark that the presence of the star 
cannot be neglected when modeling phenomenologically the observations of high- 
mass microquasars. It seems very likely that these systems cannot be understood 
as just accretion/ejection systems to first, or even zero, order of approximation, but 
require a much more complex approach including hydrodynamical and magnctohy- 
drodynamical simulations of jets and their environment, a consistent treatment of 
acceleration and emission of particles, and the physics of OB stars and their winds. 

5.3. Prospects 

The future instruments like MAGIC II, HESS II, CTA, and the already on-flight 
GLAST will allow to make a big step in our understanding of microquasars and 
gamma-ray emitting binaries. The better performance at low energies of all these 
instruments will allow to study for instance the possible development of electromag- 
netic cascades (in case of weak magnetic fields) or the secondary single-scattering 
IC component, giving information on the physical conditions of the VHE emitter 
environment. Furthermore, radiation components below 100 GeV coming from re- 
gions invisible in TeV (e.g. due to severe absorption or maximum particle energy 
i 100 GeV), could also be investigated. The good sensitivity in the whole energy 
range would allow the study of fast variability and accurate modeling, which is of 
primary importance to understand the structure of the source, and the nature of the 
particles and processes behind the gamma-ray emission. Finally, a extension of the 
operation energy range up to ~ 100 TeV (e.g. CTA) would bring the opportunity to 
study in detail the acceleration processes with extreme efficiencies, like in the case 
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of LS 5039. 
6. Summary 

Microquasars are sources in which non-thermal processes occur. Particle accelera- 
tion is taking place in these systems, although the mechanisms involved are still 
unclear (e.g. non-relativistic and relativistic shocks, magnetic turbulence, magnetic 
reconnection, etc.). The complexity of such a multiwavelength emitters, in which 
different radiative processes may be relevant in the same energy range, and emission 
reprocessing via photon-photon absorption is common, yields the study in detail of 
the underlying physics quite difficult. Nevertheless, the detection of TeV photons 
can help to constrain the fundamental properties of the emitting region, since the 
relevant timescales are short. Because of this, Cherenkov astronomy is allowing for 
the first time to really probe the accelerator/emitter in microquasars. Cygnus X-l, 
observed by MAGIC, has been found flaring during phases when photon-photon 
absorption is expected to be very severe, giving information on the emitter loca- 
tion and the stellar magnetic field. LS 5039, detected by the Cherenkov telescope 
HESS as a periodical TeV emitter, is at the moment the best laboratory to un- 
derstand particle acceleration and radiative processes in galactic compact sources. 
LS I 61 +303, detected as well by MAGIC, shows TeV emission that varies along 
the orbit. Since geometric effects would not play in this source a role for the orbital 
variability as important as in the case of LS 5039, intrinsic properties of the emitter 
should also change. The same would apply to Cygnus X-l. 

We conclude that a leptonic scenario can explain the radiation from micro- 
quasars at very high energies, although hadronic emission, energetically very de- 
manding, cannot be discarded. In case of LS 5039 and Cygnus X-l, acceleration of 
particles and emission seem to take place in the borders of the system. Due to their 
relatively hard VHE spectra, the emitter magnetic field in LS 5039 and LS I +61 303 
should be low. Due to the magnetic field produced by an OB star, gamma-rays, al- 
though likely photon-photon absorbed, can hardly trigger efficient electromagnetic 
cascades and secondary energy may be radiated mainly via the synchrotron process. 

TeV microquasars are showing us that they behave in a different way from their 
extragalactic analogs, the blazars, being very much affected by the presence of the 
primary star and the the orbital motion. 
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